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ABSTRACT 


The  hull  lift  coefficient  is  one  of  the  most  important 
and  one  of  the  least  tractable  components  of  the  dynamic 
straight-line  stability  and  maneuvering  equations.  An 
analysis  of  the  basic  flow  condition  has  been  made,  and 
similarities  to  theories  and  data  for  flat  plates  have  been 
examined.  The  hull  flow  condition  is  too  complicated  to  be 
compared  directly  to  flat  plate  flow  conditions,  but  useful 
insight  still  has  been  gained.  Available  data  for  wing 
sections  and  sets  of  systematic  hull  series  have  shown  trends 
of  lift  coefficient  versus  B/L  and  Cg,  and  have  allowed 
rudimentary  analyses  of  the  effects  of  skegs,  propellers, 
and  rudders. 
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NOMENCLATURE 


2 

aspect  ratio  = span  /area  [treating  ship  hull  as  a 
double  body,  A = (2T)72TL  = 2T/L] 

ship  or  section  maximum  beam 

block  coefficient  = displaced  volume/L*B*T 

1 2 

lift  coefficient  = lift/^pU  (characteristic  area) 

drag  or  resistance  of  water-borne  body 

ship  length,  usually  Lbp,  but  also  Lwl  for  purposes 
of  this  report 

ship  draft 

lateral  force  coefficient  on  hull  due  to  sway 

1 12- 
-(LiftAj-pU  characteristic  area  + Drag/^-pU 

characteristic  area) 

angle  of  attack,  measured  from  centerline 

2 

(single  prime)  refers  to  non-dimensionalization  by  L 
(double  prime)  refers  to  non-dimensionalization  by  LT 
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I NTRODUCT I ON 


Dynamic  stability  and  turning  characteristics  are 
important  items  in  ship  design.  In  recent  years,  the 
course-keeping  problems  of  supertankers  have  made  headlines 
which  have  made  naval  architects  more  aware  of  rudder, 
skeg,  and  hull  design  for  maneuvering.  Less  heralded  have 
been  the  studies  of  turning  that  warship,  patrol  craft, 
and  high  speed  containership  designers  have  made.  In  the 
first  two  categories,  designers  have  sought  "high  maneuver- 
ability," meaning  that  they  want  the  ships  to  "turn  on  a 
dime"  yet  still  go  straight  ahead  when  directed  to  do  so. 

The  containership  designers  have  been  concerned  with  the 
"excessive"  straight  line  stability  of  their  high-powered, 
long  and  slender  ships,  since  getting  the  ships  to  turn 
requires  enormous  control  surfaces. 

In  this  report,  I had  set  out  with  the  intention  of 
studying  the  maneuvering  problem,  specifically  the  derivative 
Yv,  of  high  speed,  long  and  slender,  transom  stern  vessels. 
Although  some  very  nice  photographs  of  flow  under  a transom 
stern  patrol  craft  were  published  by  Saunders,  very  little 
quantitative  data  was  available.  The  net  result,  in  dealing 
with  the  data  available,  is  that  much  more  study  has  been 
made  of  medium  to  very  full  hull  forms.  A few  data  points  for 
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individual  destroyer-type  ships  are  plotted,  but  the 
general  trends  and  analyses  are  derived  from  the  full-form 
ships. 
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1.  THEORY  OF  LATERAL  FORCE  ON  SHIPS 

The  generalized  equations  of  motion,  as  set  forth  in 
Reference  2,  yield  a linear  mathematical  model  for  steering 
and  maneuvering  a surface  ship: 

(X-  - m)u  + XuAu  = 0 (1) 

(Y*  - m)  v + Y v + (Y  • - mx  ) r + (Y  -mu.)r  + Y^5=0  (2) 
v v r G r i o 

(N*  - mx  )v  + N v + (N*  - I )r  + (N  - mx_u. )r  + N~6  = 0 (3) 
vo  v rz  r G l 6 

The  criterion  for  dynamic  stability  in  straight  line  motion 
follows : 

Yy(Nr-mxGu1)  - Nv(Yr-mu1)  > 0 (4) 

Thus,  for  dynamic  stability,  the  important  coefficients  are 

Yv,  Y^,  Nv,  and  Nr.  The  first,  Yv = lateral  force  coefficient, 

is  not  easy  to  calculate  by  analytic  methods.  It  is  one  of 

the  most  important  coefficients  and  requires  accuracy  to  be 

useful.  Jacobs  [17,18]  has  shown  that  N can  be  estimated  as 

v 

a subsequent  operation  to  the  calculation  of  Y^.  The  value 
for  Y^  usually  is  quite  small,  hence  simplified  calculation 
procedures  suffice.  The  value  for  N^,  though  large,  also  can 


V 


be  treated  by  simplified  procedures  to  obtain  sufficient 
accuracy. 

Attempts  have  been  made  in  the  past  to  fit  empirical 
curves  for  Yv  to  existing  data.  The  most  commonly  cited 
effort  is  Jacobs'  [17,18]  use  of  the  Jones  [19]  formula  to 
treat  bare  hulls,  then  to  add  the  effect  of  large  area  of 
deadwood  aft.  The  result  was  a formula: 


+ D " 
o 


(5) 


where  3CL"/9a  = ttA  for  the  complete  ship. 

The  formula  gives  values  that  are  "of  the  right  orders  of 

magnitude  and  indicate  correct  trends."  It  can  be  seen 

from  Figure  11  that  the  formula  always  overpredicts  the 

value  of  Y.  ' even  without  the  addition  of  D '.  Inoue  [16] 
v o 

modified  the  Jacobs  formula,  making  use  of  the  fact  that 
Dq'  is  quite  small  in  comparison  to  dCL'/da: 

2 2 

Y ' = 1.6  it  (T  /L  ) plus  a correction  multiplier  (6) 

for  trim 

or 

Yv"  = • 8 ttA  (7) 

Gerritsma  et.al.  [12]  cited  this  formula  as  being  principally 
based  on  Dollay's  [6]  low  aspect  ratio  theory  and  a number 
of  empirical  allowances.  However,  it  very  much  resembles 


the  Jones  formula  with  a correction  multiplier. 

The  Jones  formula  indeed  does  a fair  job  of  predicting 
Yv' . It  takes  no  account  of  factors  other  than  aspect 
ratio  of  the  basic  hull.  This  objection  alone  certainly  is 
not  significant,  since  correction  factors  easily  could  be 
applied.  However,  a much  more  fundamental  objection  exists. 
The  Jones  formula  was  derived  from  a sound,  theoretical 
examination  of  triangular  or  delta  low  aspect  ratio  wings. 
Jones  considered  a long,  flat,  triangular  airfoil  travelling 
point- foremost  at  a small  angle  of  attack.  The  theory 
depended  on  the  expansion  of  sections  in  a downstream 
direction.  It  is  true  that  ships  expand  from  the  bow  in  a 
downstream  direction,  but  this  is  an  expansion  in  thickness, 
not  in  section  width.  According  to  the  theory,  any  part  of 
the  airfoil  surface  having  parallel  sides  would  develop  no 
lift.  Further,  sections  aft  of  the  section  of  maximum 
width  would  not  develop  lift.  Jones  used  an  assumption, 
"corresponding  to  the  Kutta  condition,"  to  show  that  the 
shape  of  the  foil  aft  of  the  section  of  maximum  width  has 
little  or  no  influence  on  the  lift  of  the  foil.  Jones 
showed  that  the  tip  of  the  foil  need  not  be  pointed;  an 
elliptic  or  similar  plan-form  having  expanding  section- 
width  in  a downstream  direction  gives  nearly  the  same 
results.  The  reason  can  be  surmised  from  the  figures  below: 
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The  evidence  of  observations  of  flow  over  the  bottom 
of  ship  model  hulls  indicates  that  the  lift  on  ships  depends 
heavily  on  the  conditions  near  the  trailing  edge.  An  ideal 
fluid  treatment  of  a wing-shaped  body  in  unbounded,  invis- 
cid  flow  results  in  the  generation  of  a large  moment,  often 
referred  to  as  the  "Munk  moment"  (21),  but  no  lift.  Exam- 
ination of  the  result  of  including  viscosity  reveals  that 
net  lift  forces  are  created  (see  Appendix  II).  "The  total 
lift,  as  a result  of  an  inertia  distribution  and  a viscous 
distribution  along  the  ship  length,  is  generated  for  the 
greater  part  in  the  forebody,  which  means  that  the  viscous 
part  counter-balances  nearly  the  inertia  part  in  the  after- 
body. The  center  of  the  viscous  force  distribution 
therefore  lies  well  aft  of  the  center  of  gravity."  [121 
This  last  point  makes  it  quite  evident  that  the  fit  of  the 
Jones  formula  to  ship  data  is  happenstance. 

Jones  stated  that  the  theory  was  not  good  for  rectang- 
ular planforms.  He  showed  that  for  such  forms,  the  theory 
predicts  an  infinite  concentration  of  lift  at  the  leading 
edge  and  no  lift  elsewhere,  whereas  a more  accurate  theory 


would  show  some  distribution  of  the  lift  rearward.  Bollay 
introduced  a theory  that  could  treat  rectangular  planforms. 
Starting  with  a thin  flat  rectangular  plate  having  very 
low  aspect  ratio,  Bollay  treated  the  lift  problem  as  a 2-D 
phenomenon,  ignoring  leading  and  trailing  edge  effects. 

For  the  infinite  aspect  ratio  case,  Bollay  found  that  the 
wing  is  hydrodynamically  equivalent  to  a constant  distribu- 
tion of  bound  vortices  whose  axes  lie  along  the  span.  At 
the  ends  of  the  span,  these  vortices  separate,  flowing 
downstream  along  the  resultant  velocity  vector.  Inclination 
of  the  vortices  to  the  wing  was  assumed  to  be  most  important 
at  the  separation  point,  since  the  part  of  the  vortex 
nearest  the  wing  is  the  most  influential  in  inducing  down- 

wash  at  the  wing.  The  theoretical  normal  force  coefficient 
12  2 

is  CN  = N/^pU  area  = 2 sin  u;  hence,  at  u=0,  0^=0.  For 

aspect  ratios  in  the  range  0 < A <•  1,  the  normal  force 

2 

coefficient  retains  a sin  u form,  except  that  the  slope 

3CN/<htQ  becomes  increasingly  positive,  so  that  the  coeffi- 

2 

cient  takes  on  the  form  CKI  = C.sinn  + C~sin  <x.  At  A = 1 , C., 

N 1 2 N 

vs.  u is  nearly  a straight  line  up  to  the  stall  angle. 

Bollay  stated  that  the  theory  worked  for  A=  0 and 
A=°“,  and  that  the  low  aspect  ratio  formulation  worked 
relatively  well  for  0<A<  1.  The  reasons  for  the  gradual 
breakdown  are  that  three-dimensional  effects  occur  and  that 
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leading  edge  and  trailing  edge  effects  change  the  lift 
distribution.  Figures  13  and  15  show  this. 

The  aspect  ratio  range  for  ships  is  about  ,05<A<  .2, 
which  is  well  within  the  range  of  applicability  of  the 
Bollay  theory.  Unfortunately,  ships  do  not  resemble  flat 
rectangular  plates.  The  fore  portions  usually  do  have  a 
rectangular  planform,  while  the  sterns  often  are  cut  away. 

The  important  difference  is  the  thickness  variation  along 
the  length  and  the  shape  of  transverse  sections.  By  virtue 
of  several  factors,  including  blunt  bows,  widely  varying 
sterns,  beamy  sections,  large  Reynolds  number  flow,  and  bow 
wave  effects,  the  end  effects  become  significant.  Even  flat 
rectangular  plates  being  operated  at  the  free  surface  in 
water  are  not  directly  amenable  to  treatment  by  Bollay 's 
theory.  Bertrand  [3],  testing  a small  plate  with  A=  .26, 
in  fully  submerged  flow  in  a propeller  tunnel,  verified 
Bollay's  theory.  But  Gerritsma,  et.al.,  testing  a 3-meters 
long  plate  with  A=  .114  at  the  free  surface  in  a towing 
tank,  obtained  a lift  coefficient  about  half  that  predicted 
by  Bollay's  theory.  Part  of  the  discrepancy  may  be  attri- 
butable to  measurement  inaccuracy,  as  admitted  by  Gerritsma, 
et.al.  More  likely,  the  differences  in  testing  conditions 
account  for  very  different  vortex  and  boundary  layer 
characteristics  due  to  cavitation,  aeration,  and  free  surface 


effects.  Crabtree,  as  cited  by  Jacobs  [17],  noted  that  a 
flat  rectangular  plate  showed  a pronounced  suction  peak  in 
the  pressure  distribution  near  the  leading  edge,  with  a 
consequent  steep  adverse  pressure  gradient  and  laminar 
boundary  layer  separation.  Jacobs  further  noted  that  "since 
the  size  and  form  of  the  separation  region  or  'bubble'  has 
a large  effect  on  the  lift,  Tsakonas  cautioned  against  the 
use  of  the  experimental  measurements  on  flat  plates  to 
predict  forces  in  the  case  of  the  ship-wing  analogy."  (This 
caution  would  seem  to  apply  to  Jones'  theory  as  well  as  to 
Bollay's  theory.) 

The  net  result  is  that  flat  plate  theories  and  experi- 
mental results  have  to  be  modified  to  suit  the  realities  of 
flow  around  ship  hulls.  The  bow  and  stern  sections  of 
merchant  ship  hulls  (Mariner,  Series  60,  etc.)  begin  to 
resemble  flat  plates  with  sharp  or  squared-off  edges. 
Transom-stern  or  cutaway-stern  ships  such  as  destroyers  have 
less  resemblance  to  plates  at  their  sterns.  The  differences 
between  plates  and  ship  hull  mid-bodies  are  another  matter 
entirely.  However,  such  differences  should  not  be  inter- 
preted as  being  of  such  nature  as  to  make  the  flat  plate 
theories  useless.  What  is  necessary  is  some  way  to  bridge 
the  gap  between  the  plates  and  hulls;  to  show  how  the  hulls 
make  not  very  good  low  aspect  ratio,  thick  wings.  To  dismiss 
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the  flat  plate  data  entirely,  to  use  instead  a formula  with 
an  entirely  unrelated  theoretical  origin  simply  because  it 
runs  close  to  the  data  for  some  ship  hulls,  is  undesirable. 

It  is  more  satisfactory  to  start  from  an  idealized  concept, 
then  to  work  toward  the  reality. 

Bertrand's  study  hints  at  why  flat  rectangular  plate 
theory  is  not  useful  for  directly  predicting  the  lift 
characteristics  of  ship  hulls.  His  data  for  rectangular 
and  delta  plates  showed  some  decrease  in  CT  when  sharp  lead- 
ing  and  side  edges  were  rounded  off.  The  shed-vortex  sheet 
was  observed  to  break  away  from  the  bottom  edge  further 
downstream  for  the  rounded  edges  than  for  the  sharp  edges. 
This  suggests  that  the  rounding  off  more  easily  allowed  flow 
around  the  bottom  edge,  from  the  high  pressure  side  to  the 
low  pressure  side,  without  the  need  to  generate  tip  vortices 
at  the  forward  end.  This  effect  would  be  magnified  by  a 
ship  hull,  since  bilges  are  round  and  the  width  across  the 
bottom  increases  from  the  bow  to  the  parallel  mid-body. 

Saunders  [22)  presented  photographs  (Figure  15)  of  the 
underwater  portions  of  ship  hulls  to  illustrate  the  flow  at 
the  shell.  It  is  easy  to  see  that  the  rounded  bilges  and 
wide  expanses  of  bottom  plating  cause  the  flow  to  be  anything 
but  two-dimensional  across  the  span  (in  this  case,  from  free 
surface  to  keel) . There  is  considerable  cross  flow  being 
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swept  downstream  by  the  forward  velocity,  so  that  the  flow 
on  each  side  is  changed  considerably  from  the  waterline  to 
the  keel.  In  the  fore  portions,  there  is  little  evidence 
of  the  trailing  vortices  that  would  be  formed  along  the 
bottom  of  a flat  plate.  Hence,  little  viscous  lift  is 
produced  by  the  fore  portions  (supporting  the  assertion  by 
Gerritsma,  et.al.).  The  stern  skegs  or  deadwood  present  a 
very  different  flow  problem.  They  are  like  flat  plates 
themselves  and  do  much  to  help  generate  trailing  vortex 
systems.  It  is  tempting  to  observe  that  skegs  and  large 
areas  of  deadwood  aft  look  like  good  cases  for  application 
of  the  Jones  formula  for  delta  wings.  They  are  delta  shapes, 
fitted  to  the  hulls  with  the  hulls  acting  as  groundboarcs  so 
that  effective  aspect  ratios  are  double  the  measured  aspect 
ratios.  However,  it  is  not  apparent  from  the  photographs 
that  the  skegs  act  as  isolated  lifting  surfaces.  The 
• photographs  show  that  when  the  flow  has  reached  a skeg, 

substantial  changes  in  free-stream  flow  have  occurred.  The 
shed  vortex  system  has  been  started  by  the  keel.  On  hard 
chine  vessels,  the  chine  itself  can  be  expected  to  act  as  a 
bottom  edge  on  a very  low  aspect  ratio  lifting  surface,  with 
the  consequent  generation  of  its  own  vortex  systems.  The 
Jones  formula  may  be  of  some  use  for  predicting  the  lift 
generated  by  a small  skeg  on  a large,  full-bodied  ship,  such 


O 
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as  a tanker,  but  the  interactions  with  hull  flow  make  the 
Jones  formula  inappropriate  for  calculating  the  lift  due  to 
large  skegs  or  deadwood  areas  on  slender  ships  such  as 
destroyers.  [Note  that  Saunders'  photographs  are  of  ships 
in  fully  developed  turns,  being  tested  on  a rotating  arm. 
The  angle  of  attack  varies  along  the  length,  starting  from 
about  zero  at  the  bow  and  increasing  in  a downstream  direc- 
tion. However,  Y and  3CT/3a  are  defined  for  small  angles 
of  attack  for  the  whole  ship.  Thus,  the  photographs  are 
not  representative  of  the  flow  at  the  forebody  for  genera- 
ting Yv«  A uniform  angle  of  attack  should  generate  more 
cross  flow  forward  and  should  generate  the  beginnings  of  a 
shed  vortex  system  further  forward  than  shown.] 


i 
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2.  RESULTS  OF  BOLLAY'S  THEORY  FOR  LOW  ASPECT  RATIO 
FLAT  RECTANGULAR  PLATES  l 6 ] 


F = force  on  plate,  resultant  of  lift  and  drag 
N = normal  force,  component  of  lift  resolved  perpendicular 
to  plate 


N/L  = cosa 


1 . ..2, 


Cjj  * normal  force  coefficient  = N/^-p  U area 


2 

At  A = 0 , = 2 sin  a 

At  A=  .5,  is  nearly  linear  with-  a 

At  A=  1.  , CN  is  linear  with  a for  0°  < u < 20° 


= t — (2sin  a/cosa)  _ . n 

3a  3a  ' a = 0 , A = 0 

o 


3 2 

4sinu  + 2sin  a/cos  a 


for  A>0  takes  on  increasingly  positive  values, 


O (0  £0  30  40  so  o 40  £0  30  40  SO 

« ( DCURCCS)  o(  (D£GlR£.ES) 

bOLLAY  THEORY 


Figure  1:  Normal  Force  Coefficient  versus  Angle  of  Attack 
a for  Various  Aspect  RatioN  Flat  Plates 
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EXPERIMENTAL 


THEORETICAL 


(coefficient  x aspect  ratio)  plot  as  straight  lines,  while 
the  non-linear  low  aspect  rati~>  and  Dollay  theories  are  well 
separated. 


Aspect 

CN  e 3o° 

3Ct, 

CN  @ 25° 

3Ct, 

2 TT 

Ratio 

3a 

da 

1 + Va 

0 

0 

0 

0 

0 

0 

.0333 

.14 

.267 

.08 

.183 

.103 

.134 

.34 

.649 

.20 

.458 

.393 

i 

.20 

.42 

.802 

.33 

.756 

.571 

.35 

.62 

1.18 

.48 

1.10 

.936 

.50 

.90 

1.72 

.72 

1.65 

1.257 

.66 

1.05 

2.01 

1.00 

2.29 

1.559 

1.00 

1.33 

2.54 

1.28 

2.93 

2.094 

1.25 

1.40 

2.67 

1.55 

3.55 

2.417 

1.5 

1.50 

2.86 

1.60 

3.67 

2.693 

2. 

1.67 

3.19 

1.65 

3.78 

3.142 

°N 

at  30°  and 

CM  at  25° 

N 

are  measured 

from  lines  tangent 

to  the 

curves  at  a 

= 0. 

The  plots  of 

3CL/3a  versus  A have  been  non-i 

dimensional- 

ized  on 

jpUL’T  rather  than 

1 2 

2 f>  U L because  the 

first  shows 

up  the 

differences 

much  better.  The  formulae  of 

the  form 

Figure  3: 


Lift-Slope  ;)Cl"/;Li  versus  Aspect  Ratio 
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Within  the  ship  aspect  ratio  range,  Bollay's  plots 
(Figure  1)  of  CN  versus  a are  distinctly  non-linear.  This 
causes  no  problem  with  Yv'  for  the  stability  equation,  since 
Yv*  = OCl'/3«0)  + D ' which  is  defined  for  small  perturba- 
tions in  a about  ct  = 0,  in  Phase  1 and  early  Phase  2 of  a 
turn.  However,  the  non-linearity  will  affect  the  use  of  Yv ' 
in  the  linear  maneuvering  equation,  where  substantial  values 
of  a are  used  in  later  Phase  2 and  Phase  3 of  a turn  (Figure 
16)  . 

Depending  on  the  preferences  of  the  authors,  versions 
of  the  non-linear  maneuvering  equation  include  the  terms 

AY  = Y v + Y , i v I v I + ... 

V V V 1 1 


or 


AY  = Y v + Y v3  + 
V vvv 


in  the  Taylor  expansions  for  forces  [4], 

The  Bollay  data  indicate  that  the  first  expression  gives  a 
better  fit  to  flat  plate  results.  However,  a lack  of  data 
for  ship  hulls  makes  it  impossible  to  choose  one  form  or  the 
other  for  the  purpose  of  fitting  the  data  better.  The  non- 
linearities  introduced  by  hull  geometry  probably  complicate 
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the  problem  so  much  that  any  higher  order  between  2 and  3 
would  give  acceptable  results. 

Landweber  and  Johnson  120)  also  cited  a formula, 

L = C^a  + C2O? 

for  bodies  of  revolution,  but  the  exponent  could  have  been 
3 as  easily  as  it  was  2. 
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3.  WING-SHIP  ANALOGY 


For  the  wing-ship  analogy,  the  hull  is  treated  as  a 

very  low  aspect  ratio  fin  having  very  large  area.  Although 

the  hull's  performance  as  a fin  is  quite  poor,  the  area  is 

large;  hence,  force/unit  area  x area  becomes  a significant 

quantity.  As  with  all  wing  theories,  the  primary  parameter 

for  evaluating  different  lift  coefficients  for  different 

2 

hulls  is  the  aspect  ratio,  span  /area.  The  free  surface  of 
the  water  acts  as  a groundboard  or  mirror,  effectively 
doubling  the  span's  effect.  Hence,  span  = draft  x 2. 
Assuming  that  hulls  have  nearly  rectangular  planform,  espec- 
ially when  rudders  and  skegs  are  included,  area  = draft  x 2 
x length.  Then  aspect  ratio  = A = 4T^/2TL  = 2T/L. 

Other  parameters  are: 
thickness/chord  = B/L 

taper  ratio  = C /Cg  = tip  chord/root  chord  = 1 

sweepback  angle  = angle  that  trace  of  chord  midpoints  makes 

with  vertical  ~ 0 

mean  section  sh«»pc  = waterplane  at  load  waterline 
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4.  BASIC  WING  THEORIES  FOR  ELLIPTICAL  PLANFORMS 


For  infinite  aspect  ratio. 


Cg  = 2ua 


where  C^  = section  lift  coefficient. 

The  Lanchester-Prandtl  wing  theory  (1)  applies  to  wings 
having  elliptical  spanwisc  distributions  of  lift  and  aspect 
ratios  larger  that  about  3.  To  relate  finite  and  infinite 
aspect  ratios  at  constant  lift  coefficient. 


a * aQ  + (Cl/tiA) 


where  a = angle  of  attack  for  finite  A 

aQ  = angle  of  attack  on  a section,  infinite  A 
CL  = wing  lift  coefficient 


a = uo  + (V"A)  = ao{1  + A) 


In  the  linear  range  of  lift  coefficient  versus  angle  of 
attack , 


CL  " (8CL/3a)o  , OCL/3a)  = CL/a 


9 
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Ct  - OC£/3ao)ao  , OC£/3aQ)  * C£/aQ 


(12) 


At  CL  m C£ 


(3CL/3a)a  = OC£/3ao)aQ 


(13) 


(3CL/3a)ao(l  + J)  = (3CjL/dao)aQ 


(14) 


(3C  /3a)  (1  + J)  = (3C£/3aQ) 


(slopes) 


(15) 


(CL/a)  (1  + £)  - <C£/a0) 


(16' 


At  a * a. 


CL  “ C£*a^“o*  + 


,,  2.  2na 

^i/1  + A = 1 + 2/A 


(17) 


The  Jones  theory,  for  A less  than  .5,  is 


CL  “ 2<l  A 


(18) 


Botwcen  . 5<  A<  3.,  tip  vortices  give  a downward  compon- 
ent to  the  inflow  velocity,  thereby  decreasing  the  inflow 
angle  a.  The  expressions  for  finite  aspect  ratio  lift 
coefficient  are  applicable  only  to  flat  plates  with  ellipti- 


Q 
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5.  EXAMINATION  OF  WING  SECTION  PROPERTIES 


Although  the  wing-ship  analogy  treats  the  hull  as  a 
low-performance  fin,  it  is  possible  that  section  properties 
have  similar  effects  on  hulls  and  on  wings.  To  pursue  this 
possibility,  the  effects  of  maximum  section  thickness  ratio 
(y/c)  and  position  of  maximum  section  thickness  (x/c)  were 
examined.  The  results  were  compared  to  data  for  three 
series  of  ships,  and  a general  agreement  was  found  between 
the  effects  of  y/c  on  wing  sections  and  the  effects  of  B/L 
on  hulls. 

From  section  characteristics  [1,15]  (Figure  4), 

C"  = 3C"/3o  = 2tt  [1. 4 - £U.94  + . 8 £]  (19) 

L»  ii  C C 

(.03  < y/c  < .11) 

x/c  * LCB  from  bow/L  (20) 

y/c  = ~B/L  (21) 

Thus, 


3Clh/3«  = 2 n ( 1 . 4 - ^£^1  ( . 94  + . 4 |] 


(22) 


(.06  < B/L  < .22) 


No  ship  data  has  been  found  to  prove  or  disprove  the 
validity  of  the  factor  for  LCB  location.  Only  the  factor 
for  section  thickness  can  be  examined  in  this  report. 


CL*  = VjPV2L2 


C " = L/ipV2LT 


C ' = c " — 
L L L 


A model  series  (12]  specifically  designed  to  yield  data 
on  variation  of  maneuvering  characteristics  with  B/L  had  a 
slope  of  3Cl"/9c»  versus  B/L  of  about  1/2tt,  rather  than  0.4 
from  the  section  data.  Other  model  series  were  not  consist- 
ent, but  the  general  trend  was  to  support  a positive  slope 
within  each  series,  with  roughly  the  same  coefficient. 


O 
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DA  T C 0 M 


Deal gnat  ion 
NACA 

Table 
4.1.1-A 
Cia(per  *) 

C^per  rad) 

From  Abbott 
& von  Dounhot'  f 
Clipper  rad) 

y/c 

X 

0006 

.108 

6.188 

6.215 

3.001 

30 

0009 

.109 

6.245 

6.011 

4.501 

30 

63-006 

.112 

6.417 

6.215 

3.000 

35 

009 

• 111 

6.360 

6.071 

4.500 

35 

012 

.116 

6.646 

6.322 

6.000 

35 

015 

.117 

6.704 

6.322 

7.500 

35 

018 

.118 

6.761 

6.366 

9.000 

35 

021 

.118 

6.761 

6.548 

10.500 

35 

64-006 

.109 

6.245 

6.112 

2.995 

40 

009 

.110 

6.303 

6.388 

4.490 

40 

012 

.111 

6.360 

6.548 

5.981 

40 

015 

.112 

6.417 

6.388 

7.482 

35 

018 

.111 

6.360 

6.279 

8.979 

35 

021 

.110 

6.303 

6.433 

10.481 

35 

65-006 

.105 

6.016 

6.051 

2.998 

40 

009 

.107 

6.131 

6.031 

4.496 

40 

012 

.110 

6.303 

6.112 

5.997 

40 

015 

.110 

6.303 

6.279 

7.498 

40 

018 

.100 

5.730 

6.011 

8.999 

40 

021 

.112 

6.417 

6.548 

10.499 

40 

66-006 

.100 

5.730 

6.112 

3.000 

45 

009 

.103 

5.901 

6.112 

4.499 

45 

012 

.106 

6.073 

6.112 

6.000 

45 

015 

.105 

6.016 

5.821 

7.495 

45 

63A010 

.105 

6.016 

6.011 

4.995 

35 

64A010 

.110 

6.303 

6.011 

4.995 

40 

Data  for  PC /3n  Dependence  on  thickness/chord)  and  location 
of  Maximum  Section  Width 


LCNGlT  h 


Figure  4:  Lift-Slope  3C  "/3a  versus  ti/L,  LCD,  for  DATCOM  Section 
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6. 


Y versus  B/L 


Tests  by  Gerritsma,  Beukelman,  and  Glansdorp  [12]  on  a 

special  set  of  Series  60  hulls  showed  a definite  increase 

of  Y ' with  increase  of  B/L  beyond  a value  of  about  .05 

(L/B  = 20) . Even  after  the  effects  of  drag  were  subtracted, 

the  tests  still  showed  an  increase  of  3CT '/3a  with  increase 

of  B/L.  The  slope  for  Yv'  versus  B/L  was  about  .018-  .014/. 2 

= .02  at  Fn  = .15.  The  slope  for  3CT '/3a  versus  B/L  was 

Lt 

approximately  .01. 

The  MARAD  Series  [13]  showed  little  or  no  increase  of 
Y ' with  increase  of  B/L.  Examples  were: 

.00670  - .00645/. 200  - .167  = .0075  and 
.0126  - .0124/. 200  - .181818  = .011. 

The  original  Jacobs  Series  60  [18]  without  propellers 
or  rudders  actually  showed  decrease  of  Yv'  with  increase  of 
B/L.  The  slope  was  -(.01387  - . 01264) /. 1666  - .125  = -.03. 
However,  the  Jacobs  and  Eda  [11]  Series  60  for  the  same 
hulls  with  propellers  and  rudders  produced  a slope  of 
.01643  - .01509/. 1666  - .125  = .032. 


; 


[ 


\ 
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7 • Yv*  versus  Cp  or  Cp 

Some  small  increase  in  Yv'  is  to  be  expected  for 
increase  in  Cfi,  just  from  the  increase  in  drag  on  a fuller 
form  hull.  However,  the  drag  contribution  to  Y ' is  small 
anyway,  so  a small  addition  may  not  be  noticed.  Another 
contribution  to  Y ' arises  from  the  change  in  hull  shape. 

The  MARAD  Series  [13]  best  illustrates  the  effects  of 
CQ  variation  within  a consistent  series  of  hull  forms  (Figure 
6).  When  Cfi  increases,  L^/L  (parallel  midbody/total  length) 
increases  proportionally.  This  may  cause  an  increase  in 
Yv'  because  it  increases  the  length  over  which  midship 
section  bilge  radius  occurs  and  it  increases  the  length  of 
bottom  having  midship  section  bottom  width.  Fore  and  aft  of 
the  parallel  midbody,  the  curvature  from  flat  bottom  to  sides 
increases  in  radius,  and  the  width  of  flat  bottom  tapers  to 
zero  at  the  ends.  Two  effects  then  augment  each  other.  The 
greater  extent  of  minimum  bilge  radius  makes  the  hull  more 
like  a thick  foil  with  small  radius  edges.  The  greater 
extent  of  full  width  bottom  creates  more  skin  friction 
resistance  to  cross-flow,  since  such  flow  will  have  greater 
distances  to  traverse  fore  and  aft  of  the  former  ends  of 
parallel  midbody. 

The  series  of  elongated  bodies  of  revolution  [20]  also 


i 


O 
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have  increasing  lift  with  increasing  Cn , probably  because  the 

D 

fuller  shape  has  more  skin  friction  resistance  to  transverse 
flow,  and  presents  more  lateral  area  within  the  confines  of 

d x 4. 

The  results  of  several  series,  when  compared  series 
against  series,  show  that  gross  increases  in  C_  with 
accompanying  changes  in  other  characteristics,  generally 
lead  to  decreases  in  Y ' . The  full  form  tanker  hulls,  such 
^s  the  MARAD  Series,  the  HyA  tankers  [5,8],  and  the  Glansdorp 
and  Pijfers  tankers  [14],  lie  well  below  the  Eda  destroyers, 
and  below  the  Series  60  hulls  on  the  Y^'  versus  A plots. 
Furthermore,  the  finer-formed  series  are  less  consistent  in 
their  dependence  of  Y ' on  Cn.  This  may  occur  because  the 
increases  in  CB  decrease  the  slenderness  of  the  ends,  cut 
down  on  deadwood  area,  and  decrease  lift  generation  at  the 
ends,  especially  the  sterns. 
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8.  PROPELLER  EFFECTS 

Propellers  are  known  to  contribute  to  directional 
stability.  Their  effect  on  Y ' is  measurable,  as  can  be 
seen  from  comparison  of  powered  and  non-powered  model  test 
results.  One  form  has  been  suggested  from  submarine  analysis 
by  Abkowitz  (under  the  assumption  that  the  propeller  is  oper- 
ating near  the  point  of  maximum  efficiency,  9n/9J=0,  behind 
a streamlined  body) : 

(1  - "E»  <2^>  <kt- J |tT)  <26> 

d = propeller  diameter 
J = advance  coefficient  V /nd 

A 

2 4 

Kt  = thrust  coefficient  = T/pn  d 
L = ship  length  (Lbp) 
n = propeller  efficiency 
n = RPM 

p = density  of  seawater 
T = thrust  of  propeller 
V = speed  of  ship 

= speed  of  advance  of  propeller 
wf  = wake  fraction  = (V-VA)/V 

Unfortunately  for  this  study,  the  data  necessary  to 
calculate  AY^ ' are  not  available  in  the  data  sets  used. 
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Thus  the  formula  is  useful  only  as  a suggested  form  without 
experimental  verification  for  ship  hulls. 

A formula  of  this  type  would  be  desirable  for  evalua- 
ting the  differences  in  Yv'  attributable  to  various 
propeller  parameters,  such  as  power  per  propeller  and  number 
of  propellers.  It  is  desirable  to  have  a breakdown  of 
contributions  to  Y ' rather  than  to  have  some  simple  formula 
for  overall  Y ' which  is  invariant  with  changes  in  basic 
parameters  other  than  aspect  ratio. 
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Sample  Calculation  for  Propeller  in  Free  Stream 


AY  ' = 
v 


1. 8d‘ 


(1  - 


wf  > W (kt  - 


J — T) 

9J  1 


L = 450  ft. 

T = 27  ft. 

d = .7  T 

Fn  = .20 

V = .20  x /32 . 2 x 31>0  = 24  ft/sec. 

= 14.2  knots 

EHP  = 2800  hp  = I^V/550  = (1  - t)TV/550 


J = VA/nd  = 24  x .7/1  x 27  x 
Kt  = T/pn2d4  = EHP  x 550/(1  - 
= 2b00  x 550/(1 
= .3 


assume 


-.3 


7 = .9 
t ) V p n2d4 

.2)  x 24  x 2 x l2  x (27  x .7)  4 


AY 


= -1.8 


(27  x .7) 
4 5CT2 


( 1 - .3)  ( 


2 -n  x . 6 


) 1.3  - .9  x (-.3)  ) 


| 

4 


= .00034 


Q 
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The  propeller  operates  i >»  water  which  has  been  dis- 
turbed by  the  passage  of  the  hull,  and  in  general  the  water 
around  the  stern  has  a forward  motion  acquired  by  its 
interaction  with  the  movement  of  the  ship.  The  forward- 
moving  water  is  called  the  wake,  and  it  results  in  the 
propeller  advancing  relative  to  the  water  at  some  speed 
(speed  of  advance)  which  is  a lower  speed  than  the  ship 
speed  V. 

The  difference  between  ship  speed  V and  speed  of 
advance  VA  is  called  the  wake  speed.  Taylor  introduced 
"wake  fraction"  = w^: 

wf  = (V  - VA)/.V  (27) 

V=  V(1  - w,)  (2  8) 

A f 

The  propeller,  when  developing  thrust,  accelerates  the 
water  ahead  of  it,  thereby  lowering  the  pressure  around  the 
stern  and  increasing  the  velocity  of  the  flow  around  the 
stern.  The  effect  of  the  propeller  in  inducing  an  inflow 
velocity  reduces  the  forward  wake  somewhat  below  the  nominal 
wake  for  bare  hull. 

When  a hull  is  towed,  there  is  an  area  of  high  pressure 


i 


O 


over  the  stern  which  has  a resultant  forward  component 


r 
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reducing  the  total  resistance.  With  a self-propelled  hull, 
the  pressure  is  reduced  over  some  of  the  stern  area  by  the 
action  of  the  propeller  in  accelerating  the  water  flowing 
to  it.  The  result  is  that  the  forward  component  of  the 
pressure  is  reduced,  the  resistance  is  increased,  and  the 
thrust  necessary  to  propel  the  ship  at  a given  speed  is 
increased.  The  common  practice  is  to  view  this  increase 
in  resistance  as  a decrease  in  the  thrust  available  at  the 
propeller,  so  that  although  the  screw  provides  thrust  T, 
only  thrust  is  available  to  overcome  resistance.  The 
"thrust  deduction  fraction"  = t: 


t = (T  - R^/T 


Rt  = (1  - t)' 


(1  - t)  = thrust  deduction  factor 


The  net  result,  the  ratio  of  work  done  on  the  ship  to 
that  done  by  the  screw,  is  hull  efficiency  nH: 


nH  = PE/PT  = ¥/TVA  “ (1  " t)/(1  ‘ W) 


(31) 
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EFFECT  OF  BLOCK  COEFFICIENT  ON 
WAKE  FRACTION  (w)  AND  THRUST  DEDUCTION  FRACTION  (t) 


Series  60  (from  [9,25]) 


LCB 

2 . 5 A 

2.5A 

1 . 5 A 

1 . 5 A 

2.0A 

. 5A 

. 5A 

. 5A 

CB 

.60 

.65 

.60 

.65 

.70 

.60 

.65 

.70 

^'W 

.285 

.310 

.246 

.263 

.316 

.230 

.252 

.274 

t 

.179 

.206 

.177 

.167 

.204 

.154 

.155 

.161 

1.148 

1.151 

1.093 

1.130 

1.164 

1.098 

1.129 

1.156 

LCB 

. 5F 

. 4F 

. 5F 

. 5F 

. 8F 

1.4F 

1.5F 

1.5F 

1.5F 

CB 

.60 

.65 

.70 

.75 

.80 

.65 

.70 

.75 

.80 

w 

.218 

.247 

.271 

.367 

.363 

.229 

.255 

.311 

.362 

t 

.140 

.160 

.180 

.227 

.330 

. 136 

.156 

.171 

.238 

nH 

1.10 

1.116 

1.125 

1.221 

1.053 

1.121 

1.133 

1.203 

1.195 

LCB 

2. 5F 

2.6F 

2.5F 

3. 5F 

3.5F 

CB 

.70 

.75 

. 80 

.75 

.80 

W 

.259 

.297 

.355 

.288 

.355 

t 

.164 

. 170 

.200 

.176 

.170 

1.128 

1.181 

1.240 

1.157 

1.287 
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Wake  fraction  and  thrust  deduction  fraction  generally 
increase  with  C^.  The  hull  efficiency  generally  increases 
with  Cfi  since  w is  more  affected  than  t.  Hence,  the 
anticipated  increase  in  resistance  for  higher  CQ  is  some- 
what  offset  by  increased  hull  efficiency.  But  since  t 
increases  with  C^,  the  effect  of  the  propeller  in  reducing 
pressure  over  the  stern  and  pulling  the  streamlines  more 
toward  the  propeller,  increases  with  Cn . The  effect  of  the 

13 

streamline  modification  may  be  to  increase  the  lateral  area 
over  which  the  transverse  velocity  components  due  to  side- 
slip can  generate  lift,  thereby  increasing  transverse  lift 
and  the  coefficient  Y 

v 

« 

This  may  help  to  account  for  the  large  augmentation  to 
Y ' by  propeller  action,  as  found. by  Glansdorp  and  Pijfers 
for  full  tanker  models.  The  AY^'  that  they  found  was  much 
larger  than  the  AYv*  predicted  from  analysis  of  propellers 
in  a free  stream.  The  effect  of  the  thrust  deduction 
phenomenon  may  be  an  explanation  for  the  discrepancy. 
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9.  RUDDER  EFFECTS 


Jacobs  Series  60,  Model  2,1,1 


L = 5.0  ft. 


L/B  * 7.0 


B = .714  ft. 


B/T  = 2.68 


T = .267  ft. 


L/T  = 18.75 


Rudder  span  = .200  ft. 


Area  = .021  ft.' 


Rudder  chord  = .105  ft. 


Aspect  ratio  = 2span/chord = 3.81 


.2  x .105/5.0  x .267  = .0157 


For  this  value  of  A,  3CT/3a  - 2ir/(l+  2/A)  = 4.12  or  4. 

L 


Lift  = CL  x ip  AU2 


CL  = OCL/3a)a  - (3CL/3a)  (V/ug) 


AY  ' = (3CT'/3a)  + D ' 

V Li  U 


but  Dq'  is  negligible. 


Therefore,  using  AY'  = CL' 


AYv'  = CLVV 


TP  sMt-T  X T X L [U  X (1  - w)  ]2 


AY  . = (i£L)  ( V_)  (1)  Z^TxL1 
AYv  ( 3a  MuJ  V 1aT2„ 

S P L U 


= 4 x 


|p  .2  x .105  [U  x (1  - w)  ]2 


U x (1  - w) 


1 7277 

jP  1 u 


4xlx.2x.l05x.7 


= .00235 


Model 

Span 

Chord 

L 

T 

A 


1,1,1 

.2 

.105 

^ 

.267 

3.81 


2.1,1 

.2 

.105 

.267 

3.81 


2,1,2 

.2 

.167 

5 feet 
.267 
2.4 


2,1,3 

.2 

.080 

.267 

5. 
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3,1,4 

.2 

.105 

.267 

3.81 


3CL/3a 

4.1 

4.1 

3.4 

4.5 

4.1 

AV 

.0024 

.0024 

.0032 

.0020 

.0024 

Model 

4,1,1 

5,1,1 

6,1,1 

7,1,1 

8,1,1 

Span 

.2 

.164 

.258 

.2 

.2 

Chord 

.105 

.105 

.105 

.105 

.105 

L 

5 

•feet 

*■ 

T 

.267 

.2175 

.345 

.267 

.267 

A 

3.81 

3.12 

4.91 

3.81 

3.81 

3CL/3a 

4.1 

3.8 

4.5 

4.1 

4.1 

AV 

,0024 

.0018 

.0034 

.0024 

.0024 

Jacobs  * 

Models,  Simplified  Calculation  for 

Rudder 

Effect  on  Y 
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EFFECTS  OF  RUDDERS  AND  PROPELLERS 

Comparing  a series  of  bare  hulls  versus  a series  of 
hulls  with  propellers  and  rudders,  all  running  at  Fn  = .20, 


Jacobs  Models 

2,1,1 

2,1,2 

2,1,3 

3,1,1 

4,1,1 

Y ’ 

V 

.01627 

.01659 

.01563 

.01643 

.01509 

2 

2 

2 

3 

4 

Y ’ 

V 

.01307 

.01307 

.01307 

.01264 

.01387 

AY  1 

V 

.00320 

.00352 

.00256 

.00379 

.00122 

Calculated 
Rudder  AY  ' 

V 

.0024 

.0032 

.0020 

.0024 

.0024 

Calculated 
Propeller  AY^' 

.00034 

.00034 

.00034 

.00034 

.00034 

Total 

.0027 

.0035  • 

.0023 

.0027 

.0027 

Jacobs  Models 

5,1,1 

6,1,1 

7,1,1 

8,1,1 

Y ’ 

V 

.01130 

.02669 

.01787 

.01723 

5 

6 

7 

8 

Y ' 

V 

.00943 

.02172 

.01531 

.01365 

AY  ' 

V 

.00187 

.00497 

.00256 

.00358 

Cairo  I a ted 
Rudder  AY  ’ 

V 

.0018 

.0034 

.0024 

.0024 

Ca leu  1 a ted 
Propeller  AY^' 

.00034 

.00034 

.00034 

.00034 

Total 

.0021 

.0037 

.0027 

.0027 

-48- 


With  one  exception,  the  simple  rudder  calculation 
accounts  for  70%  or  more  of  the  total  AYv*  measured.  In 
one  case,  it  overpredicts.  The  free  stream  propeller 
calculation  is  inadequate  to  make  up  for  the  difference  in 
most  cases.  It  is  reasonable  to  state  that  the  formula  is 


not  applicable  to  cases  where  propellers  operate  aft  of 
blunt  sterns.  However,  there  is  not;  data  available  to 
examine  the  applicability  of  the  formula  to  transom  sterns, 
where  the  propellers  have  nearly  free  stream  operation. 


: 


, 


I 


i 

i 


D 


10.  INTERACTION  OF  PROPELLERS  WITH  RUDDERS 


Mandel  [9]  and  others  have  shown  that  the  propeller 
race  makes  rudders  more  effective  as  ship-turning  devices. 

The  same  wake-straightening  effect  that  causes  propellers 
to  develop  transverse  lift  also  increases  the  angle  of 
attack  on  a deflected  rudder  mounted  aft  of  the  propeller. 
Additionally,  the  increased  flow  velocity  within  the  race 
helps  a rudder  to  develop  more  lift.  Proof  of  the  effect 
on  the  rudder's  ship-turning  ability  can  be  seen  in  Figure  10, 
where  increased  propeller  RPM  increases  . 

The  wake-straightening  effect  of  the  propeller  works 
to  reduce  the  rudder's  dynamic  stability  effectiveness. 
With  the  rudder  fixed  at  zero  deflection,  any  straightening 
effect  on  the  inflow  due  to  the  propeller  race  reduces  the 
angle  of  attack  on  the  rudder.  The  situation  becomes 
complicated  by  the  nature  of  the  wake  straightening. 

Because  the  fluid  flow  is  a continuous  phenomenon,  the  wake 
does  not  straighten  immediately  aft  of  the  propeller.  The 
effect  is  gradual,  increasing  with  distance  aft  of  the 
propeller.  Evidence  of  the  confounding  nature  of  the 
phenomenon  comes  from  four  sources. 

Jacobs  found  on  Series  60  model  2 that  AYv’  varied 
consistently  with  rudder  chord  variation.  The  variation 
was  non-linear;  the  slope  of  AYv'  versus  chord  length 
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decreased  with  increasinq  chord  length.  Presumably,  the 
chord  length  additions  were  made  to  the  trailing  edge,  where 
the  wake  was  more  straightened,  hence  the  additions  became 
less  useful. 

Eda  found  on  destroyer  models  that  Y • decreased  with 
increasing  Froude  number,  hence  speed.  Destroyers  have 
nearly  free  stream  flow  to  the  propellers,  so  the  increase 
in  Y ' due  to  propellers  that  occurs  for  full  form  ships 
will  not  be  as  large  for  destroyers.  The  high  power 
destroyers  have  considerable  wake  straightening  due  to 
propellers,  so  the  effectiveness  of  rudders  naturally 
decreases  with  increasing  Fn,  as  the  wake  is  further 
straightened. 

Glansdorp  and  Pijfers  used  two  hulls  and  two  sets  of 

rudders,  and  found  that  increasing  rudder  chord  decreased 

Y ’ at  corresponding  RPM.  "...increasing  the  rudder  area 

not  necessarily  gives  a better  course  stability,"  because 

other  things  remaining  fixed,  a decrease  in  ' decreases 

course  stability.  Since  the  propeller  RPM  did  not  vary, 

the  contribution  to  Yy'  due  to  a propeller  operating  behind 

a full  stern  also  did  not  vary.  The  decrease  in  Y ' then 

v 

can  be  attributed  to  the  shift  of  rudder  area  aft  when  the 
chord  was  increased,  which  shift  put  the  rudder  in  more 
straightened  wake  and  decreased  its  utility  as  a course 
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stabilizmg  device.  Since  increased  chord  length  reduced  the 
aspect  ratio  of  the  rudders,  a reduction  of  rudder  effective- 
ness could  have  occurred  regardless  of  the  flow  straightening 
effect  of  the  wake. 

Surber  124]  tested  a transom  stern  patrol  craft  with 

various  spade  rudders,  and  found  that  the  tactical  diameter 

2 

of  scheme  3 (rudder  area  = 93  ft.  , A = 1.80)  was  less  than 

2 

that  of  scheme  2 (rudder  area  = 101  ft.  , A = 1.46).  The 
improved  for  the  smaller  rudder  apparently  was  due  to 
its  higher  aspect  ratio  and  its  higher  percentage  of  area 
in  the  propeller  slipstream.  The  scheme  2 rudder  had  a 
long  chord  and  short  span,  reaching  just  below  the  propeller 
shaft.  The  scheme  3 rudder  had  a shorter  chord  and  longer 
span,  reaching  well  below  the  propeller  shaft.  The 
improvement  in  occurred  despite  the  fact  that  the  scheme 
2 rudder  area  was  slightly  aft  of  the  scheme  3 rudder  area, 
which  should  have  placed  the  scheme  2 rudder  in  more 
straightened  flow  and  should  have  improved  its  Y ^ . 

It  is  evident  that  a trade-off  arises  between  Y^  and 
Yv.  If  a designer  wants  a ship  with  small  tactical  diameter, 
he  should  place  the  rudder  at  a point  well  aft  of  the 
propeller  to  take  advantage  of  the  flow-straightening.  If 
instead  he  wants  a dynamically  stable  ship,  he  needs  as 
little  flow-straightening  as  possible  affecting  the  rudder. 
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so  he  should  place  the  rudder  as  close  as  possible  to  the 
propeller,  or  he  should  try  another  approach  by  using  twin 
rudders  placed  well  to  the  siue  of  the  propeller  slipstream. 

It  is  not  sufficient  to  depend  on  variations  of  rudder 
area  alone  to  properly  effect  changes  in  Y^  or  Y . The 
designer  will  need  to  evaluate  the  flow  conditions  to  the 
rudder  and  propeller  as  well  as  the  characteristics  of  such 
devices . 

[Note:  In  the  previous  discussions  of  rudders  and 

their  interaction  with  the  hull  and  with  propellers,  the 

linearized  equations  of  motion  from  Strom-Tejsen  [23]  have 

been  the  source  for  the  coefficients  Y and  Yr.  In  the 

v 5 

non-linear  equations,  the  Taylor  expansion  would  lead  to  use 

of  cross-coupled  coefficients  such  as  Y and  Y.  since 

r v66  6vv 

both  6 (rudder)  and  v (ship)  are  involved  simultaneously  in 
a non-linear  treatment  of  maneuvering.] 
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11.  SKEGS 


Several  authors  treat  skegs  as  isolated  lifting  sur- 
faces, using  skeg  effects  as  linear  additions  to  the  basic 
hull  stability  derivatives.  It  is  to  be  expected  that 
rudders  and  short  skegs  may  act  as  such  lifting  devices. 

They  stick  abruptly  into  the  flow,  through  or  beyond  the 
boundary  layer.  Because  of  their  small  size,  they  do  not 
have  a profound  effect  on  the  boundary  layer  nor  on  the 
vortices  being  developed  by  the  hull.  However,  large, 
faired-in  skegs  or  deadwood  do  have  a profound  effect  on 
the  basic  flow  conditions.  Saunders'  photographs  show  this 
well.  Jacobs  found  that  large  areas  of  deadwood  and  fairing 
could  augment  the  basic  hull  lift  coefficient  on  Taylor 
Series  840  hulls  by  a factor  of  2.  However,  it  is  not 
likely  that  such  a convenient  factor  would  apply  to  skegs 
on  other  hulls.  Glansdorp  and  Pijfers  found  a much  smaller 
effect  on  full  tanker  models.  Destroyer  hulls  might  have 
even  a larger  figure  because  their  basic  sterns  are  so 
smoothly  cut  away  from  somewhere  aft  of  midship  to  the 
transom  stern.  Unfortunately,  no  data  has  been  found  to 
support  this  statement.  Furthermore,  no  data  has  been  found 
to  examine  the  effects  of  small  changes  in  existing  skegs  on 
Yv<  It  is  likely  that  the  skeg  effect  on  Yv  is  linear  with 
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Skeq  changes  over  a small  range,  but  no  proof  has  been 


found. 


-55- 


12.  DATA  SETS 


Elongated  Body  of  Revolution  [20] 

The  results  of  experiments  on  elongated  bodies  of 
revolution  show  how  important  bow,  stern,  and  bottom  shape 


can 

be  for  generation  of  lift 

. The  hulls  had 

circular 

sections,  no  deadwood,  blunt 

bows,  and  tapered 

sterns . 

CP  ' 

= V/(nd2/4) l 

d = 

diameter 

l = 

length 

L'  = 

= 3CL'/3a  = L/jp  U2£2 

Y ' 
V 

= L'  + D ' 
o 

V 

= unknown,  assumed  small 

* 

d/*. 

(taken  as  aspect  ratio) 

L*  (Cp  = .55) 

L'  (Cp  = 

.1 

.0054 

.0063 

.125 

.0076 

.0091 

. 143 

.0096 

.0111 

.167 

.0124 

.0141 

.2 

.0166 

.0189 

.25 

.0233 

.0267 

The  plots  of  L'  versus  A fall  well  below  the  plots  for 


ships  and  flat  plates.  Figure  11.  This  is  to  be  expected 

because  there  are  no  lift  generation  devices  at  the  stern 

or  the  bow  and  the  middlebody  presents  a shape  that  the 

fluid  flows  under  very  easily.  With  regard  to  bow  and 

stern  lift  generation,  Landweber  and  Johnson  pointed  out 

that  such  effects  cannot  be  treated  as  independent  from  the 

effects  of  the  bare  hull.  There  is  too  much  interaction 

effect  with  the  boundary  layer  to  treat  skegs  and  deadwood 

as  isolated  lifting  surfaces. 

Landweber  and  Johnson  found  that  bodies  of  revolution 

have  a relation  between  lift  and  angle  of  attack  similar  to 

that  found  by  Bollay  for  flat  rectangular  plates:  "At  small 

angles  of  attack  ct,  the  lift  of  a body  of  revolution  is 

2 

expressible  in  the  form  L = C^nt  + C2<*  . It  is  well  known 

that  a very  long  cylinder  at  an  angle  of  attack  is  subject 

2 

to  a normal  force  proportional  to  sin  a,  an  effect  which 
may  be  explained  by  assuming  that  the  longitudinal  and 
transverse  components  of  the  incident  flow  act  independently. 
[Bollay' s two-dimensional  flow  across  the  span]  It  appears 
reasonable  to  suppose  that  this  effect  accounts 
quadratic  term  at  small  angles  of  attack." 


for  the 
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Special  Series  60 

Gerritsma,  Beukelman,  and  Glansdorp  tests  on  Series  60  (12). 


(For  original 

Series 

60  (25), 

6.5  < L/B 

< 8.5) 

• 

B/L 

.25 

.181818 

. 14286 

.1 

.05 

0 

L/B 

4 

5.5 

7 

10 

20 

00 

L(m) 

3.048 

3.048 

3.048 

3.048 

3.048 

3.048 

B(m) 

.762 

.5542 

. 4354 

. 3048 

.1524 

.006 

T(m) 

.1742 

.1742 

.1742 

.1742 

.1742 

.1742 

CB 

.70 

.70 

.70 

.70 

.70 

1.0 

A = 2T/L 

.114 

.114 

.114 

.114 

.114 

.114 

-Yv* (Fn  = .15) 

.018 

.017 

.016 

.0145 

.014 

.015 

-Yv’ (Fn  = .2) 

.0185 

.0176 

.0175 

.015 

.014 

.016 

-Yv' (Fn  = .3) 

.0245 

.023 

.0207 

.0176 

.0145 

.016 

i 

Norrbin  formula:  Y ' 

= -1.69ti 

2 2 
(tvit) 

- .08  (BT/L2)Cb 

.01734 

.01734 

.01734 

.01734 

.01734 

.01734 

.00080 

.00058 

. 00046 

.00032 

.00016 

.00001 

Y ’ 

V 

.01814 

.01792 

.01780 

.01766 

.01750 

.01735 

Slope  of  Yv'(Fn=  .15)  is  about  (.018-  .014)/(.25  - .05)  = .02 
Slope  of  3Cr '/3a  is  assumed  to  be  about  .01 


W 


Figure  5:  Y ' versus  B/L  for  Special  Series  60 
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Effect  of  Drag,  D 

D = |p  V2  S Cd(R,  B/L) 

CD  = CF  + CR 
Cp  = Cp(R) 


For  Special  Series  60, 


Fn 

V = Fn/gL 
Re  = VL/v 


.15 

2.69 

2.69  x 106 
.0040 


.2 

3.59 

3.59  x 106 
.0038 


. 3 

5.38  ft/sec 
5.38  x 106 
.0036 


where  Cp  is  non-dimensionalized  on  surface  area. 

Roughly,  surface  area  = LTv. 

To  non-dimensionalize  on  projected  area,  LT,  Cp  must  be 

2 

multiplied  by  it.  Further,  to  non-dimensionalize  on  L , Cp 

r 

must  be  multiplied  by  T/L. 


Cp ' = Cpir  (T/L) 

C 1 

F 

.00072 

.00068 

.00065 

Approximate  CR* 

V 

.0001 

.0003 

.0009 

.0008 

.001 

.0016 

Both  CR’  and  surface 

area  are  increasing  functions  of  : 

Thus,  CD  is  a function  of  speed  and  B/L  (CD(R,B/L)  as  above). 
The  slope  of  CD'  versus  B/L  is  positive,  and  3CL’ /3a = Yy ' - CD' 
should  have  lesser  slope  than  Yv'. 
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Bare  Hull  from  Jacobs  Series  60  (18] 


MODEL 

2 

3 

4 

5 

6 

7 

D " 
o 

0.015 

0.017 

0.014 

0.017 

0.015 

0.019 

Y " 

6 

0.245 

0.237 

0.260 

0.217 

0.315 

0.287 

i " 

p 

0.230 

0.220 

0.246 

0.200 

0.300 

0.268 

B/T 

2.68 

3.12 

2.34 

3.28 

2.07 

2.68 

B/L 

0.1428 

0.1667 

0.125 

0.1428 

0.1428 

0.1428 

L/B 

7.0 

6.0 

8.0 

7.0 

7.0 

7.0 

L/T 

18.75 

18.75 

18.75 

23.0 

14.5 

18.75 

LCG/L 

0.515 

0.515 

0.515 

0.515 

0.515 

0.505 

CB 

0.6 

0.6 

0.6 

0.6 

0.6 

0.7 

S 

0.614 

0.614 

0.614 

0.616 

0.614 

0.713 

CA 

0.9772 

0.9772 

0.9772 

0.9740 

0.9772 

0.9818 

A 

0.1067 

0.1067 

0.1067 

0.0870 

0.1379 

0.1067 

Y ’ 

V 

.01307 

.01264 

.01387 

.00943 

.02172 

.01531 

-D  1 

o 

.00080 

.000910 

.00075 

.00074 

.00103 

.00101 

-L  1 

V 

0.01227 

0.01173 

0.01312 

0.00870 

0.02069 

0.01429 

2 

3 

4 

5 

6 

7 

8 

0.021 
0.256 
0.235 
2.68 
0. 1428 
7.0 
18.75 
0.475 
0.8 
0.807 
0.9913 
0.1067 

.01365 

.00112 

0.01253 

8 


A-2T/L  0.1067 


0. 1067 


0. 1067 


0.0870  0.1379  0.1067 


0. 1067 


000/0-  o 


ASPECT  R^ATlO 


1 


n 
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Hull  with  Propellers  and  Rudders  from  Jacobs  Scries  60  [18] 


MODEL 

1.1.1 

2,1,1 

2.1,2 

2,1.3 

3,1,1 

4.1,1 

5,1,1 

6,1,1 

7,1,1 

8,1,1 

CB 

.6 

.6 

.6 

.6 

.6 

.6 

.6 

.6 

.7 

.8 

B/L 

.1333 

.1428 

.1428 

.1428 

.1667 

.125 

.1428 

.1428 

.1428 

.1428 

A “ 2T/L 

.1067 

.1067 

.1067 

.1067 

.1067 

.1067 

.0870 

.1379 

.1067 

.1067 

L/T 

18.75 

18.75 

18.75 

18.75 

18.75 

18.75 

23.00 

14.50 

18.75 

18.75 

-D  M 
o 

.015 

.015 

.015 

.015 

.017 

.014 

.017 

.015 

.019 

.021 

-D  ' 

0 

.00080 

.00080 

.00080 

.00080 

.00091 

.00075 

.00074 

.00103 

.00101 

.00112 

-Y  " 

8 

.255 

.305 

.311 

.293 

.308 

.283 

.260 

.387 

.335 

.323 

-Y  ' 

V 

.0136 

.01627 

.01659 

.01563 

.01643 

.01509 

.01130 

.02669 

.01787 

.01723 

Rudder 

Area 

.021 

.021 

.0334 

.016 

.021 

.021 

.01722 

.02709 

.021 

.021 

Eda  Series  60  Controllabil 

ity  Models  [11] 

MODEL 

60B2 

60 

60R1 

60R2 

60B1 

60B3 

60H1 

60H2 

70 

80 

Rudder 

Area 

.021 

.021 

.033 

.016 

.021 

.021 

.017 

.027 

.021 

.021 

-Y  " 

8 

.270 

.305 

.311 

.293 

.308 

.283 

.260 

.349 

.324 

.354 

-Y  ' 

V 

.01440 

.01627 

.01659 

.01563 

.01643 

.01509 

.01130 

.02407 

.01728 

.01888 

Cg  Series 

60 

70 

80 

L/T  Series 

60H2 

60 

60H1 

L/B  Series 

60B3 

60 

60B1 

Rudder  Series 

60R2 

60 

60R1 

I 


© 


HyA  Largo  Tankor  Model s , with  Propc l lors  and  Rurido rs 


MODEL 

81000  m3  [5] 

190000  m 

B 

33.82  m 

47.5  m 

T 

12.92  m 

18.5  m 

L 

232.86  m 

305  m 

CB 

. 80 

. 84 

B/T 

2.62 

2.57 

B/L 

.145 

.156 

A 

.111 

. 121 

Y ' 


01360 


01530 
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Simulated  Self-Prope 1 led  Destroy or  |7 ] 


B = 46  ft. 

T = 14  ft. 

L = 440  ft. 

CB  = .494 
A = 4000  tons 
A = .0636 
B/L  = .1045 
B/T  = 3.3 

Yv'  = .10149 

T>P 


This  is  data  for  a fictitious  destroyer  escort  "study  ship", 
synthesized  from  realistic  data. 
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Navy  A0177  Fleet  Oiler  [10] 
(fully  appended  and  propelled) 


L = 550 
B = 88 
T = 32.2 
B/L  = .16 
B/T  = 2.733 

C = 27400  x 35/550  x 88  x 32.2  = .615 
B 

2T/L  = .1171 

Y • = -.016046 
v 


Mariner  [9,23] 

L = 528 

B = 76 

T = 27 

A = 19000 

CD  = .61 

B/L  = .144 

B/T  = 2.8 

2T/L  = A = .1023 

Y ' = .011604 
v 
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Eda  Destroyers,  Self-Propelled  (11) 


Ship 

DL-2 

S31  031  S41  S51 

L/T 

34 

A 

.059 

B/T 

3.41 

B/L 

.1004 

CB 

.51 

Rudder 

Area 

102.5  x 2 

V 

.01049  .009205  .00866 

S31 

sonar  dome , 

Fn 

* . 

20 

031 

no  dome , Fn 

= . 

20 

S41 

sonar  dome. 

Fn 

- . 

34 

S51 

sonar  dome. 

Fn 

= . 

40 

051 

no  dome , Fn 

= # 

40 

DD-692 

051  Fn  = .15  Fn  = .38 

27.4 
.073 
2.89 
.1055 
.53 

not  listed 


.00819  .007424  .008591  .009389 
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Glansdorp  and  Pijfers  - Full  Tanker  Models  [141 


L = 310  m 
B = 47.16  m 
T = 18.90  m 
V = 235000  m3 


B/L  = .152 
B/T  =2.5 
A = 2T/L  = .122 
V = 232690  m3 


Rudders  Balanced  Balanced  Mariner-type  Fixed  Portion  Twin 


I 

II 

III 

IV 

V 

ak 

2.87 

2.3 

3.36 

2.3 

3.37 

AreaR 

75.4 

115 

78.8 

115 

38x2 

freaR 

Lx  T 

.013 

.02 

.013 

.02 

.0065x2 

SHIP 

PROP  & RPM 

RUDDER 

SKEG 

-Y  ' 
V 

01 

.01562 

02 

Y85 

I 

.01938 

03 

II 

.01733 

04 

Y67 

II 

.01836 

05 

Y85 

II 

.01821 

06 

Y102 

II 

.01878 

Ml 

.01438 

M2 

Yes 

.01510 

M3 

Y85 

.01600 

M4 

Y85 

Yes 

.01717 
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SHIP 

PROP  & RPM 

RUDDER 

SKECi 

-Y  ' 

V 

M5 

III 

.01559 

M6 

III 

Yes 

.01619 

M7 

Y85 

III 

.01716 

M8 

Y85 

III 

Yes 

.01825 

M9 

IV 

Yes 

.01663 

M10 

Y85 

IV 

.01762 

Mil 

Y67 

IV 

Yes 

.01768 

M12 

Y85 

IV 

Yes 

.01813 

M13 

Y102 

IV 

Yes 

.01829 

Ml  4 

V 

Yes 

.01722 

M15 

Y85 

V 

Yes 

.01832 

Bare 

Hull 

Fully 

appended  hull, 

average 

01 

Ml 

0 

M 

with  skeg 

-Y 

.01562 

.01438 

.01900 

.01750 

.01820 

v 


Rudder  Propeller  RPM 


No  Prop 

67 

85 

102 

I 

.00072 

.00163 

.00275 

.00372 

+ 

II 

.00124 

.00250 

.00305 

.00377 

* 

III 

.00094 

.00217 

.00277 

.00330 

O 

IV 

.00127 

.00275 

.00344 

.00411 

A 

V 

.00057 

.00169 

.00208 

.00237 

0 
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Contribution  AYv'  due  to  propellers  is  .00162  which  is  about 
10  times  the  number  expected  from  Abkowitz's  formula  for 
propeller  in  a free  stream.  Apparently  there  is  much  more 
effect  of  a propeller  on  flow  ahead  of  and  behind  the  pro- 
peller plane  on  a full  form  ship  than  there  is  in  a free 
stream. 


1 


I 


f 

I ’ 


I 


■ 


% 


i 
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13.  CONCLUSION 


A disappointingly  small  amount  of  data  was  available 
for  analysis  of  transom  stern  vessels.  Many  more  data  were 
available  for  full  form  tankers,  so  most  of  the  results  are 
more  applicable  to  that  class  of  ships. 

Although 


V = -K3cL/3ot)  + Do'] 


Inoue  and  others  have  found  that  the  two  factors  follow 
similar  trends,  and  have  found  it  expedient  to  combine  the 


two.  Besides,  most  datasets  present  Y ' alone,  and  may  or 


may  not  include  resistance  data. 


Figure  11  shows  that  Jacobs'  version  of  the  Jones 


formula,  (tt/2)A  , always  overpredicts  the  value  of  Y ' , 


although  it  may  come  close  to  the  destroyer  data  if  propellers 
and  rudders  were  removed  (Figure  12).  Inoue's  more  recent 
version  runs  through  the  data  for  Series  60,  while  it  under- 
predicts for  fine  destroyers  and  overpredicts  for  full 

\ 

tankers  (Figure  12) . 

Bollay's  theory,  expectably,  well  overpredicts  the  value 


of  Y ' for  bare  hull.  From  the  beginning,  it  was  obvious 


that  ship  hulls  are  not  good  approximations  to  flat  rectangu- 
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lar  plates.  But  the  insight  into  the  basic  flow  phenomena 
of  vortex  sheet  generation  has  been  valuable,  especially  in 
light  of  Saunder's  photographs,  showing  vortex  formation 
along  a patrol  craft  underbody. 

My  conclusion  is  that  the  data  should  be  left  in  the 
graphical  form,  and  that  data  points  should  be  added  as  they 
become  available.  Eventually,  a function 

Yy'  = f (A,  B/L,  Fn,  CB,  ...) 

will  become  evident;  but  for  now,  the  graphical  form  with 
supplementary  graphs  for  B/L  and  CQ  dependence  is  clearest 
and  most  useful. 

Much  more  basic  experimental  work  is  desirable.  Such 
things  as  V form  versus  U form  bows;  systematic  variations 
of  skeg  length;  variations  of  rudder  size,  aspect  ratio, 
and  location  relative  to  propeller;  and  LCB  variation  have 
not  been  examined  because  there  are  no  data  available.  More 
flow  studies  like  those  of  Saunders  and  Bertrand  are 
necessary  for  basic  understanding  of  the  vortex  phenomena. 
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4*  Jacobs  Series  60,  Figure  11 
■f*  Eda  Series  60,  Figure  12 

X Mariner 

Glansdorp  and  Pijfers  tankers 
O MARAD  tankers 

A Gerritsma,  Beukelman,  and  Glansdorp  Special  Series  60 

0 Destroyers 

HyA  tankers 
^ Navy  fleet  oiler 


t. 


boll  ay 

rntoKY 


ASPCCT  R>\TIO 


and  Rudders 


FISH-EYE  VIEW  OF  THE  CROSS  FLOW  UNDER  AN  ESCORT-SHIP  MODEL  WHEN  TURNING 


Figure  15: 


Excerpts  from  Saunders'  Hydrodynamics  in  Ship  Dor.  ion 
Volume  til  (22 J . 
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Figure  15  (continued) 
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Figure  16: 


Phase  of  a Turning  Maneuver  [9] 


- M 


-85- 


REFERENCES 


1.  Abbott,  I.H,  and  von  Doenhoff,  A.E. , Theory  of  Wing 
Sections , Dover,  Inc.,  New  York,  1959. 

2.  Abkowitz,  M.A.,  Stability  and  Motion  Control  of  Ocean 
Vehicles , MIT  Press,  Cambridge,  Massachusetts,  1969. 

3.  Bertrand,  J.E.,  "Experimental  Investigation  of  Low 
Aspect  Ratio  Aerofoils,"  M.S.  Thesis  in  Ocean  Engineer- 
ing, MIT,  Cambridge,  Massachusetts,  January  1972. 

4.  Bishop,  R.E.D.,  Burcher,  R.K.,  Parkinson,  A.G.,  and 
Price,  W.G.,  "Oscillatory  Testing  for  the  Assessment  of 
Ship  Manoeuvrability,"  Tenth  Naval  Hydrodynamics 
Symposium,  MIT,  Cambridge,  Massachusetts,  July  1974. 

5.  Bjerregaard,  E.T.D.,  letter  to  Abkowitz,  M.A.,  from 
Skibsteknisk  Laboratorium,  Lyngby  Denmark,  July  1975. 

6.  Bollay,  W. , "A  Theory  for  Rectangular  Wings  of  Small 
Aspect  Ratio,"  Journal  of  Aeronautical  Science,  Vol.  4, 
No.  7,  May  1937. 

7.  Brown,  S.H.  and  Alvestad,  R. , "Simulation  of  Maneuvering 
Characteristics  of  a Destroyer  Study  Ship  Using  a 
Modified  Nonlinear  Model,"  Journal  of  Ship  Research, 

Vol.  19,  No.  4,  December  1975. 

8.  Chislett,  M.S.,  letter  to  Abkowitz,  M.A.,  from  Hydro-Og 
Aerodynamisk  Laboratorium  (HyA) , Lyngby  Denmark,  Septem- 
ber 1973. 

9.  Comstock,  J.P.  (editor) , Principles  of  Naval  Architec- 
ture , The  Society  of  Naval  Architects  and  Marine 
EngTnecrs,  New  York,  1967. 

10.  Cox,  G.G.  and  Motter,  L.E.,  "Prediction  of  Standard 
Maneuvering  Characteristics  for  a Naval  Auxiliary  Oiler 
(A0177  Class),"  Report  SPD-624-01,  Naval  Ship  Research 
and  Development  Center,  Bethesda,  Maryland,  June  1975. 

11.  Eda,  H. , "Maneuvering  Performance  of  High  Speed  Ships 
and  Catamarans , " Report  SIT-DL-74-1626 , Davidson 
Laboratory,  Stevens  Institute  of  Technology,  Hoboken, 

New  Jersey,  January  1974. 


-86- 


12.  Gerritsma,  J.,  Beukelman,  W.  , and  Glansdorp,  C.C.,  "The 
Effects  of  Beam  on  the  Hydrodynamic  Characteristics  of 
Ship  Hulls,"  Tenth  Naval  Hydrodynamics  Symposium,  MIT, 
Cambridge,  Massachusetts,  July  1974. 

13.  Gertler,  M.  and  Kohl,  R.E.,  "Resistance,  Propulsion, 
and  Maneuverability  Characteristics  of  MARAD  Systematic 
Series  for  Large  Full-Form  Merchant  Ships,"  Report 
7370-1,  Hydronautics , Inc.,  Laurel,  Maryland,  November 
1974. 

14.  Glansdorp,  Ir.  C.C.,  and  Pijfers,  J.G.L.,  "The  Effect 
of  Design  Modifications  on  the  Natural  Course  Stability 
of  Full  Tanker  Models,"  Report  335,  Laboratorium  Voor 
Scheepsbouwkunde , Technische  Hogeschool  Delft,  The 
Netherlands,  November  1971. 

15.  Hoak,  D.E.  (project  engineer) , USAF  Stability  and 
Control  DATCOM,  Flight  Control  Division,  Air  Force 
Flight  Dynamics  Laboratory,  Wright-Patterson  Air  Force 
Base,  Ohio,  October  1960  (revised  January  1976). 

16.  Inoue,  S.,  "The  Determination  of  Transverse  Hydrodynamic 
Nonlinear  Forces  by  Means  of  Steady  Turning,"  Eleventh 
International  Towing  Tank  Conference  Proceedings,  Tokyo 
Japan,  October  1966. 

17.  Jacobs,  W.R.,  "Method  of  Predicting  Course  Stability  and 
Turning  Qualities  of  Ships,"  Report  945,  Davidson 
Laboratory,  Stevens  Institute  of  Technology,  Hoboken, 

New  Jersey,  March  1963. 

18.  Jacobs,  W.R.,  "Estimation  of  Stability  Derivatives  and 
Indices  of  Various  Ship  Forms,  and  Comparison  with 
Experimental  Results,"  Report  1035,  Davidson  Laboratory, 
September  1964. 

19.  Jones,  R.T.,  "Properties  of  Low-Aspect-Ratio  Pointed 
Wings  at  Speeds  Below  and  Above  the  Speed  of  Sound," 
Report  835,  NACA. 

20.  Landweber,  L.  and  Johnson,  J.L.,  "Prediction  of  Dynamic 
Stability  Derivatives  of  an  Elongated  Body  of  Revolution, 
Report  C-359  (recently  declassified) , David  Taylor  Model 
Basin,  Washington,  D.C. , May  1951. 


-87- 


21.  Munk,  M.M.  , Aerodynamic  Theory  of  Airships,  W.F.  Durand, 
Editor,  Julius  Springer , l9l4 . 

22.  Saunders,  H.E.,  Hydrodynamics  in  Ship  Design,  Vol.  Ill, 
The  Society  of  Naval  Architects  and  Marine  Engineers, 
New  York,  1965. 

23.  Strom-Tejsen,  J.,  "A  Digital  Computer  Technique  for 
Prediction  of  Standard  Maneuvers  of  Surface  Ships," 
Report  2130,  David  Taylor  Model  Basin,  Washington,  D.C. 
December  1965. 

24.  Surber,  W.G.  Jr.,  "Turning  and  Maneuvering  Tests  of 
Model  4945  Representing  the  U.S.  Coast  Guard  350-foot 
High  Endurance  Cutter  with  the  "V"  Form  Bow,"  Report 
1835,  David  Taylor  Model  Basin,  Washington,  D.C.,  April 
1964. 

25.  Todd,  F.H. , "Series  60  Methodical  Experiments  with 
Models  of  Single-Screw  Merchant  Ships,"  Report  1712, 
David  Taylor  Model  Basin,  Washington,  D.C. , July  1963. 


-88- 


APPENDIX  I 

Some  Established  Characteristics  of  Hydrodynamic  Derivatives 

1)  represents  the  slope  of  the  curve  of  the  X force 
plotted  against  Au.  Direction  of  X force  is  opposite  to 
direction  of  drag  or  resistance  force.  Xu  will  be  negative 
as  long  as  the  net  drag  increases  with  speed.  For  displace- 
ment ships,  Xu  will  be  large.  For  planing  or  other  non- 
displacement ships,  at  the  point  of  planing  or  lift-off  the 
drag  can  decrease  with  an  increase  in  speed,  making  Xu  a 
positive  quantity,  briefly.  X as  a force  represents  the 
imbalance  between  drag  or  resistance  force  and  force  pro- 
vided by  propeller  thrust.  At  equilibrium  speed  uQ,  X = 0. 

2)  X^u  represents  the  force  that  the  body  experiences  in 
the  x-direction  as  a result  of  acceleration  in  the  x- 
direction.  The  inertial  reaction  force  of  the  fluid  against 
the  acceleration  that  the  body  causes  is  in  the  direction 
opposite  to  the  acceleration  of  the  body.  X^  is  negative 
and  is  5%  to  10%  of  the  mass  for  normal  ship  types. 

3)  Y^v  represents  the  linear  approximation  of  the  Y force 
resulting  from  acceleration  in  the  y-direction.  The  origin 
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of  the  force  is  similar  to  that  for  X*u.  Y*  for  both  bow 

u v 

and  stern  is  strongly  negative,  the  magnitude  being  close 
to  that  of  the  mass  for  normal  ship  types. 

4)  N^.r  represents  the  linear  approximation  of  the  moment 

resulting  from  angular  acceleration  about  the  origin  of  the 

axis  system  fixed  in  the  body.  A positive  r produces  a 

local  v = rd^  at  each  point  P(b)  forward  of  the  origin  and 

a local  v = -rd  at  each  point  P(s)  aft  of  the  origin.  The 
s 

inertial  reaction  force  of  the  fluid  against  the  acceleration 
is  a hydrodynamic  force  in  the  direction  opposite  to  the 
acceleration.  The  moments  caused  by  the  forces  add  to 
produce  a significant  negative  value  for  positive  r. 

5)  Y^r  results  from  the  imbalance  of  the  forces  which 
cause  N^.r.  But  for  Y^.r,  the  forces  on  the  bow  are  in  the 
direction  opposite  to  the  direction  of  the  forces  on  the 
stern.  Hence,  Y^.r  may  be  ( + ) , (-)  , or  in  rare  cases,  zero. 
The  magnitude  of  Y^  will  be  small.  If  the  bow  predominates, 
having  greater  pressure  distribution,  Y^  will  be  negative. 

6)  N^v  arises  from  the  same  hydrodynamic  inertial  reaction 

that  causes  Y*v.  But  since  the  bow  and  stern  forces  are  in 
v 

the  same  direction,  N^v  will  result  from  an  imbalance.  Hence 


9 
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the  magnitude  of  may  be  (+) , (-) , or  zero.  Thu  magnitude 

of  will  be  small.  If  the  bow  predominates,  will  be 
negative. 

7)  Yvv  results  from  the  steady-state  lift  that  arises 
from  the  angle  of  attack  of  the  body  as  it  moves  ahead  at 
velocity  uq  and  to  the  side  at  velocity  v.  A positive  v 
causes  forces  in  the  negative  y-direction  on  all  sections 
along  the  length  of  the  body.  These  forces  add  to  give  a 
significant  negative  value  to  Y . 

8)  Nvv  arises  from  the  same  hydrodynamic  lift  forces  that 
cause  Yvv.  But  as  with  N^v,  the  imbalance  of  bow  and  stern 
forces  causes  Nv  to  be  (+) , (-) , or  zero,  and  to  have  small 
magnitude.  If  the  bow  forces  predominate,  Nv  will  be 
negative . 

9)  Y^r  arises  from  the  lift  generated  as  the  angle  of 
attack  increases  as  the  body  moves  ahead  at  velocity  u . A 
positive  r causes  negative  Y forces  on  the  bow  and  positive 
Y forces  on  the  stern.  Y represents  the  imbalance  of  these 
forces,  hence  it  will  be  small.  If  bow  forces  predominate, 
Yf  will  be  negative. 
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10)  Nrr  has  the  same  origin  as  Yfr,  but  now  the  bow  and 
stern  forces  add  to  make  N^_  a significant  negative  value. 
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APPENDIX  II 

Munk  Moment  versus  Actual  Moment  [Empirical]  for  Submarines 

In  non-viscous  flow,  unstable  moment  for  forward  motion 
is  proportional  to  (k^  - k^)  where  k 2 and  k^  are  coefficients 
of  added  mass  in  lateral  and  longitudinal  directions. 

2 

Munk  Moment  = moment  of  instability  = (k^-k^JVpfV  /2)sin2u 

Empirical  formula  for  submarines  (roughly,  an  ellipsoid  of 
revolution)  [20] , 

Mw'  = 0.87(k2  - k^m'  = -N  J 

= M /i p L3u 
w 2 

m'  = m/i p L3  = 2V/L3  ; V = displaced  volume 


4 2 

For  ellipsoid  of  revolution,  V = jir  b a ; 2a  = L,  2b  = H 
Therefore , 
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Pressure  forces  in  potential/ 
non-viscous  flow,  due  to  2 
stagnation  points. 


In  potential  flow  with  circulation,  the  aft  stagnation  point 
moves  to  the  tail  to  satisfy  the  Kutta  condition: 


This  is  brought  about  by  circulation  flow  about  the  airfoil: 


-96- 


which,  when  combined  with  the  incoming  flow,  creates  lift, 
the  pressure  distribution  for  which  is,  typically; 


The  net  effect  is  to  redistribute  the  forces  which  cause  the 
moment  of  instability; 


and,  in  the  submarine  case,  to  reduce  the  moment  by  13%. 


Hence , 


(Lift  x moment  arm)  = (Munk  Moment)  - (Actual  Moment) 

= (Munk  Moment)  - .87 (Munk  Moment) 
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